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Abstract 

Usual modal analysis techniques are based on the Fourier transform. Due to the AT.A/ 
limitation, they perform poorly when the modal overlap n exceeds 30%. A technique based 
on a high-resolution analysis algorithm and an order-detection method is presented here, 
with the aim of filling the gap between the low- and the high-frequency domains (30% < 
// < 100%). A pseudo-impulse force is applied at points of interests of a structure and 
the response is measured at a given point. For each pair of measurements, the impulse 
response of the structure is retrieved by deconvolving the pseudo-impulse force and filtering 
the response with the result. Following conditioning treatments, the reconstructed impulse 
response is analysed in different frequency-bands. In each frequency-band, the number of 
modes is evaluated, the frequencies and damping factors are estimated, and the complex 
amplitudes are finally extracted. As examples of application, the separation of the twin 
modes of a square plate and the partial modal analyses of aluminium plates up to a modal 
overlap of 70% are presented. Results measured with this new method and those calculated 
with an improved Rayleigh method match closely. 

Key words: Modal analysis. Medium-frequency domain. Large modal overlap. Structural 
acoustics 
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1 Introduction 



In the dynamic response of a structure, three spectral domains are usually defined: 
low-, mid- and high-frequency. In general, each mode is described by a modal 
frequency, a modal damping factor, and a modal complex amplitude distribution 
(see e.g [[T]| or The low-frequency domain is characterised by distinct reso- 
nance peaks and the strong modal character of the vibratory behaviour. When the 
frequency increases, the traditional modal identification methods cannot be used: 
damping increases, resonances are thus less pronounced, modes overlap and the 
frequency-response tends to a smooth curve. In the high-frequency domain, the 
vibration can be described as a diffuse wavefield (see e.g. Il3l,[|4l or ||5l). 
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Figure 1. Scheme of individual modal resonances with the same amplitude and a modal 
overlap factor of 100% (after HI) 



The modal overlap factor /i is the ratio between the half -power modal bandwidth 
and the average modal spacing: = J ^'^^ (see e.g. H). The boundaries of 

'—^J mode 

the three spectral domains are established according to the values of ji. One could 
define the low-frequency domain as the domain of application of modal analysis 
techniques: individual modes can be distinguished. It is generally admitted that the 
modal analysis techniques based on the Fourier transform meet their limits when 
the modal overlap reaches 30% (see e.g. H or [|71); this is due to the AT.Af limi- 
tation of this signal processing method. 

It is commonly considered that high-frequency is reached for /i = 100% (see 
FiglU): the diffuse wavefield approximation becomes valid [|5l. In this spectral do- 
main, the Skudrzyk's mean- value method (flU and [[8l) identifies a structure by 
its characteristic admittance, which is equivalent to the admittance of an infinitely 
extended structure. Adding other hypotheses, it is possible to apply statistical meth- 
ods such as the Statistical Energy Analysis (SEA) H, which seeks to calculate the 
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spatial average of the response of each component of a structure by considering 
the equilibrium of power flows. Besides the diffuse wavefield of each subsystem, 
the assumptions required by SEA are that the system represents a reverberant field, 
that the input power sources are uncorrelated, and that the subsystems are weakly 
coupled ([9J and LIOJ). 



In the hope of filling the gap between the low- and the high-frequency domains 
(30% < /i < 100%), or in effect, extending the low-frequency domain, a tech- 
nique based on the high-resolution analysis algorithm ESPRIT IfTTI and the order- 
detection method ESTER lfT2ll is described here. Three examples of application 
are presented: the separation of twin modes of a square plate (local modal over- 
lap ji = 200%) and two partial modal analyses of aluminium plates up to a modal 
overlap fi = 70%. 



In this article, modal analysis is restricted to linear systems; therefore, the impulse 
response ^ (x, t) at any point located in x is expected to be a sum of complex expo- 
nentials (decaying sines): 



e(x,t)=3? 



'K/2 

J2 afc(x)e-"^*e2"J'^**+^'^^-W 

k=l 



(1) 



where K/2 is the number of modes, fk are the modal frequencies (in Hz), the 
modal damping factors (in s~^), afc(x) and <y^?fc(x) the modal amplitudes and phases 
at the point of interest. 



The free dynamics of the generalised modal displacement qk is ruled by the follow- 
ing differential equation* 



rrikqk + Ckqk + mkUijqk = (2) 



where irik is the modal mass (in kg), Ck the modal damping coefficient (in kg s ^) 
and Uk the modal angular frequency (in rad s~^). 

The modal damping factor ak (also called modal decay constant in s~^), the modal 
decay time Tk (in s), the modal loss factor rjk (dimensionless) and the modal damp- 
ing ratio Cfe (dimensionless) are related between them and to the above physical 
quantities as follows: 

ak = —= Vk = 2Cfc = Ck = 7^ (3) 

Tk 2 Uk 2mkUJk 



If A/_3dB is the same for two successive modes around /, the modal overlap /x 
becomes: 

M = xi = Af 

^ J mode ^ J mode ^ 
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In practice, the modal damping factor a and the modal local density 
estimated in average over a narrow frequency-band centered on /. 



A/ 



are 



mode 



Measured signals always contain some noise which we suppose to be additive. 
After discretisation of ^(x, t) at the sampling rate = T^^, the signal model of 
the free response of the system becomes: 



e.(x) = 



'K/2 

E 

k=l 



I...N (5) 



In order to estimate the modal parameters, a high-resolution method is applied to 
the complex signal associated to (x) . Historically, the Prony [1 3 J or the Pisarenko [il4ll 
methods rely on the resolution of a linear prediction equation. More recent tech- 
niques assume that the signal is a sum of complex exponentials added to white noise 
and project the signal onto two sub-spaces. The space spanned by a finite-length 
vector containing successive samples is decomposed into the subspace spanned by 
the sinusoids (signal subspace) and its supplementary (noise subspace). The MU- 
SIcf^lfTSH. Matrix Pencil[ll6|l, and ESPRIlFllfTTII algorithms are based on this prin- 
ciple. The latter is chosen here since it takes into account the rotational invariance 
property of the signal subspace, ensuring a more precise and robust estimation. 

In practice, the noise deviates from white noise and noise-whitening may prove 
necessary prior to analysis. A second conditioning step described by Laroche IfTTl 
consists in splitting signals into several frequency-bands: this reduces the number 
of (sub-)signal components to be estimated by ESPRIT within reasonable limits 
and is achieved by filtering the impulse response. When narrow subbands are cho- 
sen, noise-whitening usually becomes unnecessary. The next conditioning steps 
aim at reducing the length of each subband signal in order to keep the memory 
allocation low enough and the algorithm tractable in practice: each subband signal 
is frequency-shifted toward zero and down-sampled. The down-sampling factor is 
adjusted as to avoid aliasing. 

In ESPRIT, the dimensions of both subspaces must be chosen a priori and the 
quality of the estimation depends on a proper choice for these parameters. The best 
choice for the dimension of the modal subspace is the number of complex expo- 
nentials in the signal. This number is K, twice the number of decaying sinusoids. 
It is therefore advisable to estimate this number prior to the analysis. This is done 
by means of the recently published ESTER technique lfT2l . 



^ Multiple Signal Classification 

^ Estimation of Signal Parameters via Rotational Invariance Techniques 
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Figure 2. Block diagram of the high-resolution modal analysis method. 
The block diagram in (FigO describes the three main steps of the method: 

• reconstruction of the acceleration impulse response (section 12.11) : 

• signal conditioning (section 12.21) : 

• order detection and determination of modal parameters, which constitute the 
heart of the method (section [23]) . 



2 Data acquisition and signal processing 



2.1 Reconstruction of the acceleration impulse response 



A standard measuring technique in modal analysis consists in applying a pseudo- 
impulse force with an impact hammer on a structure and to measure both the ap- 
plied force and the resulting vibration, generally by means of an accelerometer. 
Taking advantage of the assumed linearity of the system, the reciprocity theorem 
is invoked in order to obtain the modal shapes: the point of excitation is varied 
while the accelerometer is kept fixed, instead of the opposite. This experimental 
procedure has been followed throughout this article. 

The analysis of free vibrations becomes a modal analysis when the response is 
normalised to the excitation of the system. The usual technique for this purpose is 
the division of the Fourier spectrum of the response by that of the excitation. In 
principle, the result is the Fourier transform of the impulse response of the system 
at the point of interest. Since our method works in the time-domain, it would be 
necessary to calculate the inverse Fourier transform of this response. In practice, 
the division of spectra proves to be dangerous for the applicability of the method: 
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quasi-zeros in the denominator introduce high- amplitude individual components in 
the ratio; they may then be transformed into quasi-sinusoids by the inverse Fourier 
transform and appear as false modal components. In our case, the normalisation 
has been achieved by reconstructing the impulse response by means of an inverse- 
filtering technique applied to the response of the system. 

The displacement g of a linear mechanical system is: 

q = gimp * / (6) 

where gimp is the impulse response. The system will be considered initially at rest 
( v(0~) = 0) in a frame of reference such that g(x, 0~) = g(x, 0"*") = at any 
position X. Without loss of generality, one may also consider gimp(O^) = 0. It 
should be noted that WimplO"*^) and v(0~^) are not zero in general. 



dt 



Denoting Laplace transforms by uppercase letters, the generic expression £ 

— /(0+) of the Laplace transform of the time-derivative of a function yields: 

The impulse acceleration response is given by: 

r = wMmp ■ F - t;(0+) 

= "'C(7imp)+t^:mp(0+)] ■ F - t;(0+) 

7 = {7imp * /} + t^imp(0+) ■ / - ^;(0+) ■ 6 (7) 



Given the measurements of the force /meas and the acceleration 7meas, the impulse 
response 7imp is estimated as follows. The first step consists in finding a finite- 
impulse-response (FIR) filter g that transforms the force signal /meas into a nor- 
malised pulse (Fig. |3]): 

/meas *g = 6p^ (8) 

Here, g stands for the impulse response of the filter in the continuous time-domain 
and 5 p represents the Dirac impulse shifted in time for causality reasons: 5 p = 



When a hammer is used to excite the system, the excitation duration is finite and 
an upper bound for the number m of samples in /meas can be given with certainty. 
After discretisation, the convolution equation Eq. ([8]) defines a system of linear 
equations. The best solution, for example in the least mean-square sense, can be 
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Figure 3. Reconstruction of one impulse acceleration response (arbitrary units), (a) Opti- 
mised filter g, with p = 276. (b) Measured force /meas ("i = 46) and reconstructed pulse 6 
shifted by p samples, (c) Measured response 7meas and reconstructed response 7imp. 



found by commonly available algorithms. We have chosen a filter with 11 m + 1 
coefficients and p = 6m. 



In a second step, Eq. (|7]) is applied to the measured values of the acceleration 
7 = 7meas(^) ^nd the force / = /meas(^)- After convolution of Eq. © by and 
substitution of /meas(i) * 5* by 5^, the result is shifted up in time by pTg. One ob- 
tains an estimation of the impulse acceleration response 7imp: 

7imp = {7meas * g}-P - ^^imp(0+) ■ S + v{0^) ■ Q-p (9) 



When the system is excited by a continuous force (no shock), f (0"*") is and the 
above expression becomes simpler. Otherwise, f (0"^) can be estimated by integrat- 
ing 7(t). In practice, it may be difficult to extract the signal from the noise in 7meas 
and obtaining a precise value of w(0+) may turn difficult. The solution consists in 
defining the origin of time slightly before the impact hammer touches the structure 
(this is generally obvious by inspection); this guarantees that ''){Ts) and t;(0"^) are 
truly 0. 



The process of retrieving the acceleration impulse response is illustrated in Fig. [3l 
The first sample of the impulse response cannot be retrieved since fimp(O^) is not 
known. If necessary, it could be reconstituted at the end of the modal analysis and 
the corresponding correction be applied to the modal amplitudes and phases. 



2.2 Signal conditioning 



2.2.1 Reduction of the number of points 

The number of operations in the ESPRIT algorithm is 0{N ^) and the computing 
duration is excessively long for a large number N of samples. Numerical insta- 
bilities may also appear. In order to overcome these problems, we adopt the pro- 
cedure proposed by Laroche [fTTl and reviewed in the introduction: band-filtering, 
frequency-shifting, and decimating. A few minor transformations are introduced. 

It is advisable to evaluate roughly the spectral density of mode^. This helps to 
define frequency-bands that contain less than say K = 2h complex components 
[[TSl . A band-pass filter between f\ and /h is designed for each band. Although not 
as efficient as IIR filters, FIR filters are preferred because their transfer function has 
no pole and therefore, does not introduce spurious modes into the signal. Various 
techniques for synthesising the filter are available. We have chosen the Blackman 
window. 

The signal is then filtered as follows. An impulse response encounters a large vari- 
ation at t = and decreases afterwards. In order to minimise the effect of the 
transient response of the filter, the signal is time-reversed prior to convolution with 
the FIR h of the filter. This does not alter its spectrum. Convolution adds a number 
of samples equal to the length of h, at the end of the reversed signal. These points 
must be removed from the beginning of the signal after it is time-reversed again 
(see below). Once filtered, only K /2 modes are kept. However this number is still 
to be determined with precision. The amplitudes and phases of the modes at the 
measured point are altered by the filtering and their transformed values are written 
hk and (pk- 

The Hilbert transform of the filtered signal is computed in order to eliminate the 
negative-frequency content of the spectrum which would cause aliasing problems 
in the next steps of the procedure. We have used the hilbert function proposed 
by Matlab®. The procedure does not include any spectrum division; the Gibbs phe- 
nomenon (very rapid oscillations) associated to the Fourier truncation done in this 
procedure is limited to the very beginning and to the very end of the transformed 
signal. Because of a very fast decay rate, it never proved problematic in practice 
(in other words: no pseudo-poles were added by the Fourier truncation). The signal 
now contains K /2 complex exponentials whose frequencies are between /i and /h. 

This signal is multiplied by exp(— 27rj /d i Ts), with i = 1, . . . ,N . This operation 
shifts the spectrum by /d which is chosen slightly less than f\. The spectrum of the 
result is now limited by // = f\ — fi and fl = /h — /d. As a matter of preference, we 

^ This may be done by mechanical reasoning or by extrapolating the low-frequency anal- 
ysis, for example. 
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have then taken the real part of this complex signal. This produces a symmetrical 
spectrum with K components between —f(^ and +f(^. 

According to the sampling theorem, the signal may now be down-sampled at a 
sampling frequency lower than F^, reducing the number of points to analyse. In 
principle, the decimating factor d could be chosen up to Fj2f(^; in practice, a safety 
margin is kept and the decimating factor that we have used was approximately 
Fs/6f(^. Requirements on the minimum number of points in the signal add other 
constraints on the decimating factor (see below). 

After decimation, time-reversing, and the removal of extra points (see above), the 
signal takes the form: 

~ N 

Si = Xi + f3i = Y.bkzl + pi i = l ... N = — (10) 

where Xi is the modal signal (to be determined), Zk = e~"* ^+'^'^jfk d j^-^ 
called poles (fk = fk — /d, ctk = Ofc), bk = dkC^'^^ are the complex amplitudes, and 
K is the number of complex exponentials to be found. ESPRIT requires that the 
number N of signal points be more than 2K. 



2.2.2 Noise whitening 

In principle, the results of the ESPRIT analysis correspond to the complex frequen- 
cies of the signal only if the additive noise (3 is white. In practical cases, the noise 
is white to first order in any narrow band, hence the interest of subband filtering 
presented above. For wide frequency-bands, including a noise-whitening step in 
the signal-conditioning procedure may improve the precision of the modal results. 
A method proposed by Badeau lfTSl consists in estimating the power spectral den- 
sity of the noise for each frequency-band and to deduce from it the corresponding 
whitening filter. The Fourier spectrum is computed first and a rankfilte^h used in 
order to smooth the spectrum. Then, the estimator of the autocovariance function 
is found by calculating the inverse Fourier transform of this filtered spectrum. A 
linear prediction on this estimator gives the coefficients of the whitening filter that 
can, finally, be applied to the original signal. 

This noise- whitening treatment did not prove necessary in the applications pre- 
sented here. 



^ In a rank filter, the data are sorted by ascending orders. The output value is the lowest 
data value, where r is the rank order of the filter. 
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2.3 Determination of modal parameters 



2.3.1 Order detection 



As mentioned above, the best choice for the dimension of the modal subspace to 
be given to the ESPRIT algorithm is K. Obviously, a larger value may also be 
chosen: some of the effective noise will be partly projected on to modal subspace, 
producing very weak or highly attenuated components. A choice smaller than K 
for the dimension of the modal subspace would introduce errors in the estimation 
of the modal components. 

In order to estimate the number of complex exponentials (that is: twice the number 
of modes) in the signal, we have used the ESTER (ESTimation ERror) procedure 
by Badeau lfT2]| which is schematically presented here. One notes that the first steps 
of this procedure are common with those of the ESPRIT algorithm IfTTll . 



The N signal data and the modal signal samples Xi are written in the form of 
Hankel matrices: 



^ Si S2 ... Si^ 



S 



. S 



N 



( 



X 



Xl X2 



\ 



X-n Xn~\-l 



X 



(11) 



nJ 



with / = N — n + 1, n being the sum of the dimensions of the signal and noise 
subspaces. 

It has been shown (see e.g. ifTTl or ifTSl ) that: 

• the estimation is optimal when n = N /3 or n = 2N/3, 

• the estimation quality is rapidly degrading outside this interval, 

• the estimation is only slightly degraded for n E [N/3, 2N/3]. 

In consequence, we have systematically chosen n = N /2. 

The correlation matrices are formed (computed in the case of i?ss): 



-Rss — ~r S Rxy. — T X X^ 

I I 



(12) 



For additive white noise with variance a^: 

E [i?ss] = i?xx + cr^I (13) 
which shows that the eigenvectors of i?xx are among those of R^s in the limit of 
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perfect estimation. 



The algorithm ESPRIT needs the K eigenvectors of i?xx to determine the poles 
{zf:}. It is now shown how to find both K and those eigenvectors. 

The eigenvalues (m = 1, ...,n) and the corresponding eigenvectors 
{wi, . . . , Wn} of i?ss are computed. It can be shown [|TT]| that 

• the eigenvalues are real and positive, 

• eigenvalues associated with the noise subspace are equal to cr^ (nearly equal for 
a non-white noise). 

Ordering the eigenvalues in decreasing order naturally selects the ones associated 
with the modal signal: in principle, K is the number of eigenvalues that verify 
Am > cr^ (see Eq. [131). The ESTER criterion presented below is more robust than 
this condition for the determination of K. 

W{p) is defined as the matrix formed by columns Wk' {wi, . . . ,Wp} with p < 
n. The matrix W^{p) is defined by removing the first line of W{p) and Wi{p) is 
defined by removing the last line of W{p). The following matrix $ and quantity 
E[p) are formed: 

^p) = W^{p)^W^{p) E{p) = W^{p)-W^{p)^p) (14) 
where is the pseudo-inverse of W . 

The ESTER criterion defines K as the highest p maximising J{p) = l/\\E{p)\\'^. 
In other words, K is found as the highest number such that E{K) approaches 0, 
which corresponds to the so-called rotation invariance of W. 

The case of a synthesised signal with 3 sinusoids and added noise (signal to noise 
ratio SNR = 50 dB) is shown in Fig. |4] (see Table [U for the modal components 
parameters). A threshold Jt is chosen (here: 10^), in correspondence with SNR 
and K is considered as the highest value of p for which J(p) > Jt (here: p = 6). 
This criterion proves to be very robust. 



Parameters of the test signal Parameters estimated by ESPRIT 



fk 


[Hz] 


2078.10 


2082.30 


2087.10 


2078.11 


2082.31 


2087.12 




[s-^] 


28.00 


31.00 


27.00 


27.96 


30.72 


27.02 


ak 




1.00 


0.80 


0.40 


1.00 


0.77 


0.40 


^k 


[rad] 




~(«-1.05) 


-0.52) 


1.56 


-1.05 


-0.54 



Comparison between true and estimated parameters of a synthetic signal. 
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.'.-I K = 6 



15 20 25 30 35 40 

Model order p 



Figure 4. Application of the ESTER criterion to a signal made of three damped sinusoids 
and additive noise (SNR = 50 dB, see Table [T] for the other parameters). The detection 
thi^eshold for ESTER is chosen to 100. The value K = 6 (corresponding to 3 modes) is 
clearly detected. 



2.3.2 Modal frequencies, modal damping factors, and complex amplitudes 

Once the number of modes K /2 has been estimated, the K first columns of W{n) 
are extracted to form W = W{K), the matrix of the eigenvectors of R^x- The 
purpose of the ESPRIT procedure is to derive the so-called poles Zk from this in- 
formation on the modal signal. The main steps are schematically recalled here (for 
a demonstration, see ifTTI ): 



The Vandermonde matrix and the diagonal matrix D are formed with the z^: 



1 


1 . 


. 1 




Z2 ■ 


■ 




^2 ■ 


Z^~ 


n-1 
1 


n-1 
Z2 


Z^~ 



\ 



Their rank is K and they verify: 



D 



Zl 



(0) 

V (0) ) 



= V^D 



(15) 



where the matrice (repectively V^) are formed by eliminating the first row 
(respectively the last row) of V"". 

The rank of W is also K and therefore, a base-change matrix C can be defined 
by: 

V'' = W C (16) 
Shifting this equation up and down yields V^' = C and = V^" 
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Using Eq.fTSlyields: 

W^^ = W^C D 



Wl W^ = C D 



(17) 



This equation, denoting a so-called rotation-invariance property of W, shows 
that the poles are the eigenvalues of Wj^ W-^. 



The frequencies and damping factors of the response signal are: 



otk 



d 



In \zk\ 



(18) 



The final step consists in the determination of the amplitudes and phases of the 
modal components. To this end, the N x K Vandermonde matrix V'^ is formed. 
The complex amplitudes bk are the best solution, in the least-mean-square sense, of 
the equation: 



bi 






bK_ 




■5 TV 



The amplitudes and phases of the response are: 

\bk\ 
\H{h)\ 
(fk = arg(6fc) - arg[iJ(/fc)] 



(19) 



(20) 
(21) 



Table \T\ shows the estimated results for the synthetic signal described above. The 
error is generally less than 1% (4% for the phase of the third component and am- 
plitude of the second component). 



3 Applications 

Partial modal analyses are shown in three cases: 

• a square aluminium plate (A) with localised damping: twin modes with yU ~ 
200%, 

• a rectangular aluminium plate (B) in the mid-frequency domain (30% < < 
50%), 

• a rectangular aluminium plate (C) in the mid-frequency domain (45% < /i < 
70%). 
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3. 1 Experiments 



A pseudo-impulse force is applied by means of an impact hammer {P. C.B. Piezotron- 
ics 086D80). The acceleration is measured with an accelerometer {Briiel & Kjcer - 
ENDEVCO, Isotron 2250A-10). In all cases, boundary conditions are kept as close 
as possible to "free-free". The point of excitation is varied (see section [2TI) whereas 
the vibration measurements are made at a single point, in the vicinity of a corner of 
the plate. Under the chosen boundary conditions, this location is not on any of the 
nodal lines. 

The signal analysis described in the previous sections is applied independently to 
each pair of measurements { /meas, 7meas } • The frequency and the damping factor of 
each mode is taken as the weighted mean of all the estimated values. Weights are 
the estimations of the amplitude at each point: this gives less importance to the less 
precise estimations in the nodal regions. 

The masses of the plates (A), (B), and (C) are respectively 0.48 kg, 5.5 kg, and 
22.5 kg. Despite its relatively low mass (0.4 grams), the accelerometer causes a 
slightly negative shift of the modal frequencies. This phenomenon was evaluated 
quantitatively on plate (A) by placing a second accelerometer with the same mass 
just opposite to the first one. A frequency drift of — 0.7 Hz was observed for the 
(2,l)-mode and of — 0.3 Hz for the (l,2)-mode, both at approximately 180 Hz (see 
section [33] ). To first order, the mass loading effect of one accelerometer can be 
corrected by adding the measured drift to the modal frequencies measured in the 
situation with one accelerometer only. For plate (B) (5.5 kg), a negative drift of less 
than 0. 1 Hz was observed for the modes of interest, around 600 Hz. For the heaviest 
plate (C) (22.5 kg) a negative drift of about 0.1 Hz was observed around 1600 Hz, 
close to the uncertainty of our method for these high frequencies (see section [331) . 

3.2 Theoretical modal determination 

Only approximate solutions are known for the frequencies and the shapes of the 
conservative modes of a thin isotropic rectangular plate with free-free boundary 
conditions. Warburton [[T9ll combined a Rayleigh method with characteristic beam 
functions to obtain a simple approximate expression. In this approach, plate modes 
are assumed to be the product of beam functions: 

W^m,n){x, y) = Xm{x)Yn{y) (22) 

where x (resp. y) corresponds to the shorter length (resp. longer) of the plate and 
Xm (resp. Yn) is the m (resp. n)-th normal mode of a beam with the same boundary 
conditions as the corresponding edges of the plate. The frequency accuracy is ex- 
cellent for plates with constrained edges but it is less so when one or more edges are 
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left free. Kim & Dickinson 11201 provide an improved approximate expression by 
using the Rayleigh method in connection with the minimum potential energy the- 
orem: the deflection Wm,nix, y) includes three terms (see [Appendix A[ ). For com- 
parison with experiments, we have retained this method since the errors on modal 
frequencies are known to be less than 1% [21] with tractable frequency expressions. 
In our experiments, the uncertainties and approximations are such (see below) that 
more precise methods (Rayleigh-Ritz method, superposition, exact series solutions, 
finite element analysis, see Hurlebaus ETI for an exhaustive comparison) are not 
necessary. 



3.3 Separation of the twin modes of a square plate (A): low -frequency and high 
or low modal overlap. 



An aluminium square plate (AU4G, 300 x 300 x 1.9 mm^) is suspended by rubber 
bands. A block of foam is glued in the centre in order to increase damping. In 
principle, modes (2, 1) and (1,2) have the same modal frequency (twin modes), 
and their modal shapes are similar under a 90° rotation (Fig. [5]). In practice, modal 
frequencies and dampings are slightly different due to imperfections in symmetry 
and isotropy. Here, the modal frequencies of the two modes are ~ 178 Hz and 
^181 Hz, corresponding to a local (apparent) modal density of ~ 3 modes Hz^^. 
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Figure 5. Twin modes of a square plate (L = 300 mm): (a) (2,l)-mode; (b) (l,2)-mode. 
Dotted lines: nodal lines. Dashed lines denote where modal analyses are performed. 

The analysis is done along one side {y = 0) and along one diagonal as shown in 
Fig. [51 Plate vibrations are damped by means of a block of foam glued in the centre. 
The modal damping factors a are ~ 20 s~^, corresponding to a modal overlap of 
~ 200%. The ESTER procedure reveals two modes in the 170-200 Hz frequency- 
band, as shown in Fig. [6]^ a). They are undistinguishable in a typical Fourier spec- 
trum (Fig. [Mb)). The estimations of the modal parameters are given in Table [2] 



With the sign of the modal phase attributed to the amplitude, the modal "signed 
amplitudes" along one side are displayed in Fig. ([7]-a) together with the theoretical 
modal amplitudes for the conservative case (dashed line). Here and in what follows, 
the measured modal shapes are normalised to a maximum of 1. The amplitudes of 
theoretical modal shapes are adjusted to yield a best fit (in the least- mean-square 
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Model order p Frequency [Hz] 



Figure 6. Analysis of the first twin modes of square plate (A), (a) ESTER criterion on the 
response signal in point 9: the value K = 4 (two modes) is detected, (b) Amplitude of the 
Fourier spectrum at the same point. The length of the "useful" signal is f» 2 s (it is masked 
by noise afterwards) and increasing amounts of zero-padding were tried: beyond a total 
length of 20 s for the analysis window (as retained here), the spectrum does not change 
appreciably, o marks: modes estimated by ESPRIT. 







(2,l)-mode 


(l,2)-mode 


Plate with extra 


/(Hz) 


177.8 


181.0 


damping 


a (s-^) 


17.4 


22.4 


Plate without 


/(Hz) 


178.1 


181.4 


extra damping 


a (s-^) 


2.6 


3.9 



Table 2 

Plate (A), with and without artificial extra damping: estimations of the modal parameters 
of the twin modes (2, 1) and (1, 2). 

sense) to the experimental data. The modal phases are given in Fig. dTJ-b). The 
modes can be considered as clearly and adequately separated in this case of very 
high local modal overlap. The differences between measured and theoretical ampli- 
tude curves of the (2,l)-mode (particularly noticeable for x > L/2) are due to the 
mass of the accelerometer placed at x = y = 0. The light mass (0.4 grams) slightly 
modifies the modal shapes. We observed that adding one similar accelerometer at 
{x = L,y = 0) removes the asymmetry of the measured modal shape. 

Without the block of foam, damping factors are around 3 s^^, corresponding to an 
overlap of ~ 30%. The estimations of modal parameters are given in Table |2l The 
"signed-amplitudes" along one diagonal are represented in Fig. [8l 

3.4 Partial modal analysis of a rectangular plate (B): mid-frequency and moder- 
ate modal overlap (30% < /i < 50%). 

The plate {AU4G, 590 x 637 x 5.2 mm^) shown in Fig.|9]is supported by four blocks 
of foam around the centre in order to ensure high damping; boundary conditions 
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Figure 7. Separation of the twin modes along one side of square plate (A), (a) Normalised 
"signed" amplitudes; (b) Phase. • mai^ks: measured (2,l)-mode. Solid line: theoretical con- 
servative (2,l)-mode. o marks: measured (l,2)-mode. Dashed line: theoretical conservative 
(l,2)-mode. 




L/10 L/5 3L/10 2L/5 L/2 3L/5 7L/10 4L/5 9L/10 L 

Position (projection on the x-axis) 

Figure 8. Separation of the twin modes along one diagonal of plate: normalised "signed" 
amplitudes. • marks: measured (2,l)-mode. Solid line: theoretical conservative (2,l)-mode. 
o marks: measured (l,2)-mode. Dashed line: theoretical conservative (l,2)-mode. 

can still be considered as "free-free". 






Figure 9. Plate (B) with the line x = where modal analysis is performed. 
The measurements are made at 33 regularly spaced points along the long side 
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(x = 0). The sampling frequency is 50 kHz. The considered frequency-band is 
520-660 Hz; the modal overlap is about 40%. In this mid-frequency region, a typ- 
ical Fourier spectrum (Fig. [TOtb)) does not exhibit well-separated modes. The re- 
sult of the ESTER procedure is shown in Fig. [13 a), revealing four modes in this 
frequency-band. 




Model order p Frequency [Hz] 



Figure 10. Plate (B): modal analysis along the long side between 520 and 660 Hz. (a) ES- 
TER criterion on one response signal (point 32): the value K = 8 (four modes) is detected; 
(b) Amplitude of the Fourier spectrum at the same point. The length of the "useful" signal 
is w 1.7 s (it is masked by noise afterwards) and increasing amounts of zero-padding were 
tried: beyond a total length of 17 s for the analysis window (as retained here), the spectrum 
does not change appreciably, o marks: modes estimated by ESPRIT. 



The modal shapes are represented by the "signed amplitudes" in Fig. \TT\ Mass 
loading creates no visible asymmetry in modal shapes of plate B (its mass is 1.4- 10^ 
times that of the accelerometer) and the negative shift of modal frequencies is about 
0.1 Hz. With help of the theoretical analysis (three-term Rayleigh method), the 
measured modes can be identified as the (3,3)-, the (2,4)-, the (4,2)-, and (0,5)- 
modes, ranking 17 to 20 in the mode series. The estimations of the modal param- 
eters are given in Table [3] together with the corresponding approximate theoretical 
modal frequencies for Young's modulus E=7.4-10^° Pa, density p=2790 kg m~^, 
and Poisson's ratio z/=0.33 as given by the manufacturer. A detailed discussion on 
the determination of theoretical modal frequencies, their dependency on material 
properties and plate geometry, and their comparison with experimental values is 
presented in the next section for plate (C). 





(3,3)-mode 


(2,4)-mode 


(4,2)-mode 


(0,5)-mode 


/ray (Hz) 


523.7 


542.7 


587.2 


645.7 


/(Hz) 


526.0 


545.8 


603.5 


640.7 


a (s-^) 


27.4 


37.6 


72.3 


52.0 



Table 3 

Plate (B): estimations of the modal parameters between 520 and 660 Hz. Top line: conser- 
vative plate treated by the improved Rayleigh method. Bottom line: experimental. 
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Figure 11. Plate (B): modal analysis along one side (normalised "signed" amplitudes). 
• marks: measured modes. Solid lines: theoretical conservative modes, (a) (3,3)-mode; 
(b) (2,4)-mode; (c) (4,2)-mode; (d) (0,5)-mode. 

3.5 Partial modal analyses of a rectangular plate (C): mid-frequency and high 
modal overlap (45% < /i < 70%). 



In order to perform modal analysis on high-order modes (~ 200) near the acoustical 
coincidence frequency, a larger plate was considered {AU4G, 1000 x 1619 x 5 mm^) 
Modes are analysed on a 10 x 10 mesh with a 1 cm grid- step. Modal frequencies 
and damping factors are determined as the weighted means of the 100 correspond- 
ing estimations. 

Two experimental setups were developed in order to ensure free-free boundary con- 
ditions for this 22.5 kg plate: suspension by six thick rubber bands glued along one 
side of the plate and suspension by two nylon lines passing trough small holes near 
the top plate edge. Both are presented in order to illustrate the sensitivity of the 
method. The experimental values of the modal frequencies are estimated with an 
accuracy of ~ 0.1 Hz (see Fig. [14] and Table [6l) in two narrow frequency-bands 
(around 1700 Hz and 2100 Hz). Theoretical values are determined as follows. 

In the frame of the Kirchhoff-Love plate theory ll22l . the modal angular frequencies 
u;m,n are given by: 

^m,n = B k^^^ (23) 

Eh^ D 

where B = ; = — - and h is the thickness of the plate. The wavenum- 

12p(l-z/2) ph ^ 

bers km,n are determined by the plate dimensions a and b and by the boundary 
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conditions. 



Since the physical parameters of the plate are not readily available with the desir- 
able precision, we have estimated the B{E, p, u, h) factor by comparing the 18"' to 
25"^ measured modal frequencies to those given by finite element simulation^ (Ta- 
ble|4l). These particular modes are chosen because they are well- separated and the 
free-free boundary conditions are well-ensured. Minimising the average of the ab- 
solute values of the relative frequency differences between experiments and FEM 
simulations yields -Bfem = 61.0 m"^ s~^. With this estimated value0, the aver- 
age relative difference in this frequency-band is ~ 0.47%. Since the finite-element 
method introduces some spurious stiffness in the simulated system, the numerical 
value of i?FEM is certainly slightly overestimated. 



/meas (Hz) 


125.9 


139.3 


141.1 


147.5 


150.6 


154.7 


160.5 


171.0 


/fem (Hz) 


125.5 


139.2 


142.2 


149.7 


151.2 


154.7 


159.5 


170.9 


/fEM /meas r\r\ 

7 X 100 

J meas 


0.30 


0.12 


0.80 


1.49 


0.38 


0.01 


0.62 


0.08 



Table 4 

Plate (C): comparison between eight modal frequencies estimated by ESPRIT (/meas) and 
calculated by a finite-element method (/fem)- 



The modal frequencies and modal shapes of the high-order modes in the two frequency- 
bands of interest (around 1700 and 2100 Hz) are calculated with the approximate 
three-term Rayleigh method, using the values estimated above for the physical pa- 
rameters. According to Reference [|2T|. the systematic error for the first modal fre- 
quencies calculated by this method is positive and less than -\-\%. 

The result of the ESTER procedure for the first frequency-band (1685-1697 Hz, 
/i ^ 45%) is shown in Fig. W2\ a) and the corresponding partial modal analysis 
results are given in Fig. [13] and Table [5l Results for the three modes detected in 
this frequency-band are reported for both suspension arrangements. Also shown in 
Fig. [13] and Table [5] are the theoretical modal shapes and modal frequencies for 
modes (10,1 1), (9,13), and (12,4) which are the 199"^ to 201^' modes. 

Measured and calculated modes match closely. The positions of the nodal lines are 
correct for the three modes. The measured modal shapes are almost identical for the 



^ Simulations are carried out with 8-node thin-shell elements (as in |[2TI ). A mesh of 
70 X 100 elements is used, coiTcsponding to « 35 points per wavelength at 200 Hz in 
the x-direction (respectively k, 30 point in the y-direction). 

^ The values provided by the manufacturer for the duraluminium properties are 
E=7.4-10^*' Pa, /3=2790 kg m~^, and z^=0.33. With these values, the value retained for 
-Bfem coiTcsponds ioh = 4.96 mm. 
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Figure 12. ESTER criterion in two frequency-bands (Plate (C)). (a) 1685-1697 Hz frequen- 
cy-band, point 1: the value K = 6 (thiee modes) is detected, (b) 2065-2110 Hz frequen- 
cy-band, point 11: the value K = 10 (five modes) is detected. 





(10,ll)-mode 


(9,13)-mode 


(12,4)-mode 


/ray (Hz) 


1695.5 


1697.0 


1703.4 


Suspension 


rub. nyl. 


rub. nyl. 


rub. nyl. 


/(Hz) 


1689.9 1690.0 


1693.1 1693.3 


1695.7 1696.6 


a(s-') 


4.2 4.6 


4.3 4.7 


3.0 3.9 



Table 5 

Estimations of modal parameters between 1685 and 1697 Hz for two suspension condi- 
tions. Top line: conservative plate treated by the improved Rayleigh method. Bottom line: 
experimental. 

two experimental setups. These results, together with the estimation of uncertainties 
(see below) illustrate the precision and reproducibility of the method. 

The values of the calculated modal frequencies (Table[5]) are systematically slightly 
larger than the measured ones by 0.2-0.5%. Since this is also the case in the 2100 Hz 
frequency-band (see below), there must be a systematic error for wich we propose 
the following explanations, (a) The value of Bfem used for the calculation of the 
modal frequencies is overestimated, (b) The improved Rayleigh method oversti- 
mates modal frequencieS (c) In Kirchhoff-Love plate theory, the rotary inertia 
and the shear effects are ignored; for the plate considered here, the correction given 
by the more precise Mindlin theory ll23ll in the uj{k) curve is around —0.5% at 
1700 Hz. 



Uncertainties reported in Table [6] for the two suspension conditions are evaluated 



According to II19II : By the Rayleigh principle, if a suitable waveform W is assumed, 
satisfying approximately the boundary conditions, the resulting frequency value is very 
near, but higher than, the true value, because the assumption of an incorrect waveform is 
equivalent to the introduction of constraints to the system. 
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Figure 13. Plate (C): partial modal analysis between 1685 and 1697 Hz, with 
fi ~ 45 %. (a), (b), (c) Theoretical modal shapes obtained by the improved Rayleigh 
method, (a'), (b'), (c') Measured modal shapes with the rubber-bands suspension, 
(a"), (b"), (c") Measured modal shapes with the nylon-lines suspension. 



according to Eq. (1241) : 



Af ^A'.?'^' 



/ / — 

with / the weighted mean of the estimated modal frequencies and Nb the num- 
ber of estimations (100 in our case). This uncertainty estimate is pessimistic since 
deviations are not weighted here. . . 



During measurements with the first experimental setup, we noticed a drift in the 
estimation of the frequencies and possibly also in the estimation of the damping 
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(10,ll)-mode 


(9,13)-mode 


(12,4)-mode 


Rubber-bands 


A/,„b 

/rub 


7.1 ■ 10-5 


9.0-10-5 


9.4 - 10-5 


suspension 


AOi-ub 
Qrub 


8.5- 10-2 


6.7-10-2 


5.3 - 10-2 


Nylon-lines 


/nyl 


4.0 ■ 10-5 


5.1 - 10-5 


2.2 - 10-5 


suspension 




8.5- 10-2 


4.8 - 10-2 


3.9 - 10-2 



Table 6 

Uncertainties on frequencies and damping factors for the three modes between 1685 and 
1697 Hz under two suspension conditions. 

factors (see Fig. [T4Ta) for the chronological representation of these estimations). 
The second suspension setup (nylon lines in small holes) appears to be more stable 
(Figs. [T4)) . The overall -1-0.4 Hz frequency-drift in the rubber-band case is larger 
than the uncertainty in the estimation of the modal frequency. The interpretation 
for the sign of the drift on frequency, for the fact that f^^b < Ayb and for similar 
observations on the damping factors goes as follows. Rubber bands add a mass to 



1696,4 
£ 1696.2 
I 1696 
S" 1695.8 

1695.6 

1695. 

1695.2, 



(a) 



(b) 



o 

^ 4.5- 

E 3,5 • ■„■ 



(c) 



C 4 

Q. 

i3.5 



1695,8 1696 169e 

Frequency [Hz] 



Figure 14. Comparison of the two suspending setups for the (12,4)-mode of plate (C). 
• marks: rubber-bands suspension, -i- marks: nylon-lines suspension, (a) Modal frequencies 
measured chronologically; (b) Damping factors measured chronologically; (c) Scattering 
of modal frequencies and damping factors. 

the system. However, rubber bands slip slightly and the added mass decreases in 
time. This is also consistent with the very slight negative drift in the damping factor 
ttrub- The ctnyi > Onjb obscrvation is interpreted by the fact that the vibrations of the 
plate are more strongly transmitted to the suspension frame by the nylon lines than 
by rubber bands. 



In the second frequency-band (2065-2 110 Hz) where modal analysis was performed. 
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the modal overlap factor is ^ 70%. Compared with the 1685-1697 Hz frequency- 
band, the important increase in damping factor (from ^ 4 s~^ to ~ 15 s~^) and 
thus in modal overlap is due to the sudden increase in acoustical radiation when the 
frequency approaches the coincidence frequency f^. For this isotropic plate, /c is 
given (see [|24ll for example) by: 



/c 



E 



(25) 



where Ca is the speed of sound in air (^ 343 ms~^). Above this frequency, the 
wavelength of flexural waves in the plate is larger than the wavelength of acoustical 
waves in air and an infinite plate radiates sound; for a finite plate, the increase in 
radiation efficiency is gradual when the frequency approaches /c (see Fig. [T6l). In 
our case, the coincidence frequency is about 2.4 kHz. 



Despite the high modal overlap factor, the ESTER procedure clearly detects the cor- 
rect number of modes (FigfT2lb)). The modal analysis results are given in Fig. [T5] 
and Table |7] together with the results of calculations for the 243'^'' to 247"^ modes, 
corresponding to the (3,21), (5,20), (13,6), (9,16), and (12,10) modal shapes. 
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Figure 15. Plate (C): partial modal analysis between 2065 and 2110 Hz, with /^t ~ 70%. 
(a), (b), (c), (d), (e) Theoretical modal shapes obtained by the improved Rayleigh method, 
(a'), (b'), (c'), (d'), (e') Measured modal shapes. 



Matching is excellent for frequency values and correct for modal shapes. The higher 
values of calculated modal frequencies (Table|7]) can be explained as in the 1700 Hz 
frequency-band. In the (12-10)-mode case, the analysed region is essentially nodal; 
the signal to noise ratio is ^ 30 dB and the method clearly meets its limits. 

Experimental results for the damping factors are displayed as a function of fre- 
quency in Fig. [TH Theoretical results in a different configuration are available for 



24 



(3,21)-mode (5,20)-mode (13,6)-mode (9,16)-mode (12,10)-mode 



/ray (Hz) 2077.3 2086.6 2096.7 2100.3 2112.6 

/(Hz) 2069.1 2075.7 2081.5 2092.1 2097.5 

a(s-i) 10.3 10.6 15.6 15.8 18.8 

Table 7 

Estimations of modal parameters between 2065 and 2110 Hz for the rubber band suspen- 
sion. Top line: conservative plate treated by the improved Rayleigh method. Bottom line: 
experimental. 



the sake of an approximate comparison: in the case of simply supported baffled 
plate and under the assumption of the diffuse wavefield, Maidanik ll25l gives an an- 
alytical expression of the average damping due to radiation. The other contribution 
to damping of an aluminium plate is due to thermoelastic losses [|26l . The damping 
model established by Chaigne et al. [|26ll gives tttherm < 0.14 s~^ for this aluminium 
plate. This value is very small compared with radiation damping in this frequency 
range of interest. It has therefore not been taken into account by the solid-curve in 
Fig. [in The main physical difference between experimental and theoretical con- 
ditions lies in the acoustical short-circuit between the front and the back of the 
plate, reducing radiation efficiency and decreasing damping factors. This is com- 
patible with the discrepancy between the measured points and the curve given by 
Maidanik. 




Frequency [Hz] 

Figure 16. Damping factors due to the radiation of an aluminium plate. Boundary con- 
ditions and radiation conditions are not the same for the experimental and for the the- 
oretical determinations, o marks: measured damping factors in the two frequency-bands 
where modal analysis was performed. Solid line: damping curve due to acoustical radia- 
tion of a baffled simply-supported plate (after ll25l ). Dashed line: (asymptotic) damping 
factor a of an infinite plate above the coincidence frequency /c (Ooo = PaCa/{ph), where 
Pa = 1.2 kg m^'^ is the density of air). 
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4 Conclusion 



The modal analysis method presented in this article resolves cases in which the 
Fourier transform meets its limits. Partial modal analyses of vibrating plates with 
high modal overlap (up to 70%) match theoretical modal predictions. This method 
may contribute to filling the gap between the low-frequency and the high-frequency 
domains where Fourier modal analysis and statistical methods respectively apply. 
The ESTER technique appears as a good tool for estimating the modal density, 
an essential parameter for the study of vibrating structures in the mid- or high- 
frequency domains. 

At frequencies larger than those presented here the results were not as satisfactory; 
this is mainly due to the signal-to-noise ratio limitation of the signal processing 
method. Moreover, the spatial resolution of the method becomes also a limiting 
factor: the uncertainty in the position of the impact-excitation 0.5 cm) becomes 
significant compared with the grid-step (1 cm). However, a larger grid-step would 
not be acceptable at the considered wavelength (about 15 cm at 2.1 kHz for plate 
(C). 

The SNR limitation can be partly overcome by the use of a continuous excita- 
tion with a signal that allows the impulse response reconstruction by deconvolution 
techniques {swept- sine technique as in [|27l , for example). 



Appendix A The three-term Rayleigh method 



According to classical plate theory (see e.g. [28J), the maximum strain energy, or 
potential energy of bending V , of an isotropic rectangular thin plate is given by: 



1 /-a 

^max TT ^ / / 

2 Jo Jo 



(A.l) 

Efv" 

where W is the modal shape and D is — ; —. The maximum kinetic energy T 

^ 12(1 -u^) 

of the plate is: 



T^.. = '-^ / W'dydx (A.2) 
/ Jo Jo 



The Rayleigh principle yields the modal angular frequency u: 



2 2K. 



= T. (A.3) 

ra rb ^ ^ 

ph / / W^dydx 
Jo Jo 
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Kim & Dickinson ll20l extend the Rayleigh method [fT9ll by considering three terms 
mW: 



W(^^n){x,y) = X^{x)Yn{y) - cX„/(x)F„(y) - dX^{x)Yn'{y) (A.4) 

where Xm{x) (respectively Yn{y)) is the m (resp. n)-th order normal modal shape 
of a beam with the same boundary conditions as the corresponding edges of the 
plate; Xm'{x), Yn'{y) are the next higher beam modal shapes, and c and d are 
constant quantities given below. In our case, boundary conditions are free-free: 
m' = m + 2 and n' = n + 2. The modal deflection Xm of a free-free beam are 
given in Eq. (1 1) of reference |[T9ll . 

By substituting Eq. (IA.4I ) into Eq. (IA.3I) and (lA.ll ), the modal angular frequency is: 

D TT^ Cjn,n + C^C*m,n+2 + d'^Cm+2,n ~ SC-E^.n. ~ '^dE^.^m + 2cdF 

pha^h^ 1 + c2 + rf2 

(A.5) 



with 



Cm,n = Gl- + G^^^ + 2{iyR^K^ + (1 - z/)J„J„ 



n^2 



= vHm{Kn + L„) + 2(1 - I^) J™M„ 
i^m,n, = -y{K^Kn + + 2(1 - u)M^Mn 

The values of Gm, H^, Jm., Km, Lm, Mm are given in references fT9| (Table 1) and 
[|20]| (Table 1). Finally, c and d can be determined by using the minimum potential 
energy theorem: 



dV dV 

'-' ' max * max 



dc dd 



d 



Gm+2,n-Fm,n -^n,m-^ 

Cm,n+2Cm+2,n ~ F"^ 
Gm,n+2En,m Em^nF 

Cm,n+2Cm+2,n ~ F"^ 
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